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Abstract 

Four  algorithms  for  simulating  flicker  FM  phase  noise  (f~3  spectrum)  are  given,  two  old  and 
two  new.  Their  Allan  deviation  and  mean  square  time  interval  error  (MSTIE)  are  examined. 
The  MSTIE  shows  that  one  of  the  old  algorithms  has  deficient  long-term  phase  deviations  on 
average.  The  Allan  deviation  does  not  reveal  this  deficiency. 


1  INTRODUCTION 

By  a  flicker  FM  model  we  mean  a  stochastic  process  that  has  a  spectral  density  (in  a  sense  that 
can  be  made  precise)  that  is  asymptotically  equal  to  const  //3  as  /  — ►  0.  Many  quartz  oscillators 
exhibit  flicker  FM  phase  noise  over  wide  intervals  of  Fourier  frequency,  the  evidence  being  Allan 
deviation  plots  that  are  approximately  flat  over  two  or  more  decades  of  averaging  time.  Therefore, 
simulations  of  systems  containing  quartz  oscillators  need  to  include  flicker  FM  generators.  Following 
are  three  classes  of  existing  flicker  FM  generation  algorithms.  They  all  have  running  time  of  order 
A  log  IV,  where  N  is  the  number  of  points  to  be  generated. 


•  Filter  cascade  generators  [1].  White  noise  is  applied  to  a  set  of  first-order  digital  filters  in 
series;  the  filters  are  so  designed  that  the  output  is  a  stationary  process  with  approximate 
spectrum  /_1  over  a  frequency  range  whose  low  end  depends  on  the  number  of  filter  stages. 
One  can  obtain  /~3  noise  by  a  cumulative  sum  operation  on  the  /_1  noise.  Because  these 
algorithms  generate  the  output  sequentially,  they  take  little  memory. 

•  Discrete  spectrum  (DS)  generators.  Complex-valued  Hernritian  white  noise  (the  Fourier 
transform  of  real  white  noise)  is  generated  in  the  discrete  frequency  domain,  multiplied  by 
/_1/2  or  /~3/2,  and  transformed  back  to  the  time  domain.  This  is  a  special  case  of  a  general 
method  for  generating  colored  noise;  the  author  has  no  citation  for  its  actual  use  to  generate 
flicker  FM,  but  Ref.  3  of  [2]  cites  a  suggestion  for  its  use. 

1This  work  was  carried  out  at  the  Jet  Propulsion  Laboratory,  California  Institute  of  Technology,  under  a  contract 
with  the  National  Aeronautics  and  Space  Administration. 
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•  Impulse  response  (IR)  generators,  discrete-time  analogs  of  Riemann-Liouville  fractional  inte¬ 
gration  (see  [2]  for  continuous-time  power-law  noise  models).  White  noise  is  convolved  with 
a  causal  filter  that  represents  fractional  summation  of  order  1/2  or  3/2.  The  Kasdin- Walter 
algorithm  [3,4]  discrete  convolution  in  N\ogN  time  by  using  the  FFT. 

Our  principal  aim  here  is  to  compare  a  DS  generator,  an  IR  generator,  and  two  new  FFT-based 
flicker  FM  generators  that  use  the  recent  method  of  circulant  embedding  for  exact  simulation  of 
stationary  processes  [5-8] .  The  circulant  embedding  algorithm  is  given  below;  one  may  think  of  it  as 
the  discrete  spectrum  method  with  a  twist.  Two  statistical  properties  are  used  for  the  comparison: 
1)  Allan  deviation;  2)  a  form  of  mean  square  time  interval  error  (MSTIE)  with  an  initial  phase  and 
frequency  removed.  Lest  this  work  be  merely  a  discussion  of  models  and  algorithms,  the  comparison 
also  includes  the  phase  residuals  of  two  precision  quartz  oscillators  that  were  chosen  for  flatness 
of  their  Allan  deviations.  The  MSTIE  (but  not  Allan  deviation)  exposes  a  deficiency  of  the  IR 
generator,  the  same  deficiency  that  a  filter  cascade  generator  has  if  it  is  not  properly  initialized  [9] : 
the  generator’s  long-term  phase  excursions  are  too  small  on  average.  Nevertheless,  as  Schmidt  also 
found  [10],  one  can  avoid  this  deficiency  by  generating  twice  as  many  points  as  needed  and  using 
only  the  second  half. 

We  will  be  working  with  discrete-time  phase  models  xn,  — oo  <  n  <  oo,  with  sample  period  1, 
normalized  so  that  their  two-sided  spectral  densities  Sx(f),  \  f\  <  1/2,  are  asymptotic  to  |27r/|-3 
as  /  — >  0.  To  simulate  samples  x  (nro)  of  flicker  FM  phase  (time)  x  ( t )  with  one-sided  frequency 
spectrum  Sy  (/)  =  /r_i/_1  and  Allan  deviation  \Jh-  \  In  4,  multiply  xn  by yj nh-iTo. 

2  TWO  FLICKER  FM  PHASE  MODELS 


Before  defining  the  four  generators,  we  define  two  flicker  FM  models,  against  which  the  generators 
can  be  compared.  Although  each  model  for  phase  xn  is  a  nonstationary  stochastic  process,  its 
second  increment  A2xn  =  xn  —  2xn_i  +  xn-2  is  a  stationary,  mean-zero,  Gaussian  process.  The 
definition  of  xn  is  ambiguous  in  that  any  constant  phase  and  frequency  can  be  added  to  it;  to  pin 
down  xn  exactly,  we  can  fix  two  values  xa  and  xy  The  two  models  differ  mainly  in  their  spectral 
densities  near  the  Nyquist  frequency  /  =  1/2  (see  Figure  1). 

2.1  FD(3/2)  (Fractional  Difference)  Model 

For  each  value  of  the  real  parameter  5,  there  is  a  process  called  FD(d)  [11,12]  with  spectral  density 
|2sin7r /p2<5  (in  a  sense  to  be  explained  for  5  =  3/2).  This  family  has  the  convenient  property  that 
if  xn  is  an  FD(<5),  then  Axn  is  an  FD(<5  —  1).  (The  frequency  response  of  the  difference  operator  A 
is  |2sin7r/|.)  In  particular,  FD(— 1/2)  is  defined  as  a  stationary,  mean-zero,  Gaussian  process  zn 
with  spectral  density  |2sin7r/|.  For  its  autocovariance  (ACV)  sequence,  sZjTl  =  E ZjZj+n,  we  have 

Sz,n=  [  e*27r/n  (2  Sin  7 xf)df=—jT^ - rr.  (1) 

Jo  7T  u  -n2) 


482 


34th  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting 


By  definition,  xn  is  an  FD( 3/2)  process  if  A2xn  is  an  FD(— 1/2)  process.  Then  xn  is  a  process 
with  stationary  second  increments,  and  Sx(f )  =  |2sin7r/|-3  in  the  following  sense:  if  H  is  a 
finite  moving-average  filter  that  contains  A2  as  a  factor,  then  Hxn  is  stationary  and  Shx  (/)  = 

|  H  (e~i2nf)  |  Sx  (/),  where  H  (z)  is  the  z-transform  of  the  impulse  response. 

It  can  be  proved  that  an  FD(3/2)  process  xn  can  be  represented  in  the  following  way: 


/  n  \  n 

xn  -  H - H - x-ni 

\  ni)  n\ 


Eh 

3= 1 


n-juj  + 


E 

j=-oo 


hn—j 


n 


h—j  -\-  h—ni_j 


n  i 


n  >  1, 


(2) 


where  ni  is  a  positive  integer,  is  a  standard  white  noise  sequence  (independent  Gaussians  with 
mean  0  and  variance  1),  and  hn  is  defined  by  the  power  series  (1  —  z)~ 3//2  =  '%2™=ohnzn",  thus, 
hn  =  (— l)n  (_3^2),  and  hn  =  0  for  n  <  0  by  convention.  One  can  show  that  hn  ~  2^/n/ir  as 
n  — >  oo. 


2.2  Sampled  PPL  (Pure  Power  Law)  Model 

Starting  with  a  continuous-time  process  x  (t)  with  stationary  second  increments  and  spectral  density 
|2-7r /|-3  for  all  nonzero  real  /,  we  sample  it  at  the  integers  to  get  a  discrete-time  process  x  (n).  Its 
second  increment,  z  (n)  =  A2x  (n),  is  stationary  and  has  ACV 

sz  (n)  =  sx  (n  +  2)  -  4sx  (n  +  1)  +  Qsx  (n)  -  4 sx  (n  -  1)  +  (n  -  2) ,  (3) 

where  sx  (t)  is  the  generalized  ACV  [13,14]  of  x  (t): 

sx  (t)  =  ^-t2ln  \t\  fort  t^O,  sx  (0)  =  0.  (4) 

Z7T 

The  spectral  density  of  the  sampled  process  is  not  \2irf\~3,  but 

OO 

£  127 T(f  +  k)\-3,  |/|  <1/2  (5) 

k=— oo 

(see  Figure  2). 


3  TWO  GENERAL  ALGORITHMS 

Three  of  the  generators  under  discussion  can  be  quickly  given  in  terms  of  two  algorithms  of  wider 
utility;  the  second  algorithm  uses  the  first.  With  small  changes,  these  descriptions  follow  Percival 
and  Walden  [5],  We  use  a  standard  complex  discrete  Fourier  transform  (called  the  FFT  from  now 
on)  and  its  inverse,  defined  for  M -point  sequences  by 
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M—l  M—l 

Ak  =  Y2  an  exp  (—i2Trkn/M) ,  an  =  —  ^  Ak  exp  (i2nkn/M)  , 

n= 0  k= 0 

respectively. 

3.1  Discrete  Spectrum  Algorithm 

Purpose:  Generate  values  of  a  real  stationary  Gaussian  process  with  a  desired  discrete  spectrum. 

Inputs:  Nonnegative  numbers  So,  . . .  ,  Sn,  where  A  is  a  power  of  2  and  Sk  is  the  desired  two-sided 
spectral  density  at  frequency  fk  =  k/  (21V). 

Outputs:  Gaussian  random  variables  zo, . . .  ,zn  such  that  E zn  =  0, 

1  N 

E zmzn  =  -  ^  5,fc,  exp  (i2vr/fc  (n  -  m)) . 

V  k=l—N 

Procedure: 

Generate  Uq,  Ui,  ...  ,  Un,  Li,  . . .  ,  hjv-i  as  independent  standard  Gaussians. 

Let  Zq  =  VSoUo,  Zn  =  sJ~SnUn- 

Let  Zk  =  y/Sk/ 2  (Uk  +  iVk),  Z2N-k  =  Z |  for  k  =  1  to  N  -  1. 

Let  the  real-valued  sequence  zq,  . . .  ,z2n~i  be  \/2N  times  the  inverse  2A-point  FFT  of  the 
Hermitian  sequence  Zq,  . . .  ,  Z2n-\- 

Keep  the  values  zo, . . .  ,  zn  (or  any  N  +  1  consecutive  values). 

3.2  Circulant  Embedding  Algorithm 

Purpose:  Generate  values  of  a  real  stationary  Gaussian  process  with  a  given  autocovariance  (ACV). 

Inputs:  Real  numbers  so>  •  •  •  ,  sn >  the  desired  ACV  up  to  lag  N,  where  A  is  a  power  of  2. 

Outputs:  Gaussian  random  variables  zq,  ■  ■  ■  ,zn  such  that  E zn  =  0,  E zmzn  =  sn-m  for  0  <  m  < 
n  <  A. 

Procedure: 

Let  sn  =  sn  for  n  =  0  to  A,  s2N-n  =  sn  for  n  =  1  to  A  —  1  (extension  of  sn  by  reflection). 
Remark:  “Circulant”  refers  to  the  covariance  matrix  that  corresponds  to  sn. 

Let  the  real- valued  sequence  5q,  . . .  ,  S2n- i  be  the  2A-point  FFT  of  sq.  . . .  ,  S2jv-i- 
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If  any  Sk  <  0,  the  method  fails.  (This  means  that  the  extended  circular  sequence  is  not 
positive-definite.) 

Use  So, . . .  ,  Sn  in  the  discrete  spectrum  algorithm  to  generate  the  zn. 


4  FOUR  FLICKER  FM  GENERATORS 

All  these  generators  produce  approximately  N  phase  values  xn .  where  A  is  a  power  of  2.  The  first 
two  are  approximate;  they  do  not  simulate  either  target  model  exactly.  The  last  two  give  exact 
simulations  of  the  two  target  models. 


4.1  DS  —  Discrete  Spectrum 

Let  fk  =  k/  (2 N).  Run  the  discrete  spectrum  algorithm  with  input  So  =  0,  Sk  =  {2n fk)~3, 
k  =  1, . . .  ,  N,  and  output  xo,  ■  ■  ■  ,  x at.  (One  could  also  set  S^  =  (2  sin  it  fk)3  to  approximate  the 
FD(3/2)  model  more  closely.) 

4.2  IR  —  Impulse  Response 

This  generator  is  an  approximate  simulation  method  for  the  FD(3/2)  model.  Its  output  is  given 
by  the  formula 

n 

xn  =  ^2  hn-jUj ,  1  <  n  <  N,  (6) 

3= 1 

where  hn  is  defined  after  (2),  and  u\, . . .  ,  un  are  independent  standard  Gaussians.  The  convolution 
is  carried  out  by  zero-padding  the  sequences  to  length  2 N,  Fourier  transforming  them,  multiplying 
the  transformed  sequences,  and  inverse  transforming  the  result.  For  details,  see  [3],  Observe  that 
(6)  is  just  one  part  of  (2). 

4.3  FD  —  Fractional  Difference 

This  generator  is  an  exact  simulation  of  N  +  3  values  xn  of  the  FD(3/2)  model.  Run  the  circulant 
embedding  algorithm  using  the  ACV  (1)  for  the  input  so, . . .  ,  sat,  and  zo,  ■  ■  ■  ,zn  as  output.  It  can 
be  proved  that  the  algorithm  succeeds.  (The  ACV  satisfies  Craigmile’s  criterion  [8].)  This  produces 
an  exact  realization  of  N  +  1  values  of  FD(— 1/2).  Then  perform  two  cumulative  summations: 

Un  =  Un—l  +  Zn- 1  for  n  =  1  to  N  +  1, 
xn  =  xn-i  +  yn- 1  for  n  =  1  to  N  +  2. 

The  initial  values  yo  and  xq  are  arbitrary,  and  may  be  set  to  zero. 
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4.4  PPL  —  Pure  Power  Law 

This  generator  is  an  exact  simulation  of  N  +  3  values  of  the  sampled  PPL  model.  It  is  identical 
to  the  FD  generator  just  described  except  that  the  ACV  (3)  is  used  in  place  of  (1).  Again,  it  can 
be  proved  that  the  algorithm  succeeds.  There  is  one  complication:  to  avoid  catastrophic  roundoff 
error  in  (3),  use  the  asymptotic  approximation 

s*(n)~~i(1  +  ^  +  i)  (7) 

in  place  of  (3)  whenever  n  >  35. 


5  COMPARISONS 

We  compare  the  four  flicker  FM  generators  with  each  other  and  with  the  phase  residuals  of  two 
quartz  oscillators  (Oscilloquartz  and  C-MAC)  that  were  compared  once  per  second  against  hydrogen 
masers.  The  test  runs  were  chosen  for  flatness  of  Allan  deviation  between  1  and  1,000  seconds2. 

Figure  1  shows  a  sample  output  of  the  four  generators  with  N  =  1024.  Also  shown  are  the  first 
1,025  phase  residuals  of  the  quartz  oscillators,  scaled  up  as  explained  in  the  next  section.  The  exact 
generators  are  both  initialized  so  that  xo  =  x\  =  0.  The  IR  output  starts  with  a  small  nonzero 
value  of  x\.  The  DS  output,  which  realizes  a  stationary  process,  has  no  special  initial  value. 

Figure  2  shows  the  spectral  density  of  the  two  target  models  along  with  the  discrete  spectrum  of 
the  DS  generator  for  N  =  32.  The  spectral  densities,  multiplied  by  (27 r/)3,  are  plotted  on  a  linear 
scale  against  /.  As  the  next  section  shows,  the  rise  in  the  sampled  PPL  spectrum  (5)  near  the 
Nyquist  frequency  is  just  right  to  make  its  Allan  deviation  exactly  flat  for  all  integral  r.  The  DS 
spectrum  actually  has  too  little  high-frequency  power  for  this  purpose,  the  FD  spectrum  too  much. 
For  many  purposes  these  high-frequency  deviations  are  of  little  concern. 


5.1  Allan  Deviation 

Figure  3  shows  the  theoretical  Allan  deviation  (lines)  and  the  square  root  of  measured  Allan  variance 
(small  dots),  averaged  over  10,000  trials,  for  the  four  generators  with  N  =  1,024.  Also  shown  are 
the  Allan  deviations  of  the  quartz  oscillators  (symbols)  as  estimated  from  their  test  runs,  which 
lasted  about  45,000  s  for  the  Oscilloquartz,  129,000  s  for  the  C-MAC.  The  oscillator  results  were 
normalized  so  that  ay  (64  s)  =  -\J (In  4)  /it  =  0.664,  the  theoretical  value  assumed  by  the  PPL  model 
for  all  t.  The  a  axis  has  an  expanded  linear  scale  to  bring  out  the  differences  among  the  plots.  The 
IR  generation  was  actually  performed  with  N  =  2, 048;  IR1  refers  to  the  first  half  of  the  generated 
sequence,  which  is  equivalent  to  IR  with  N  =  1,024;  IR2  refers  to  the  second  half. 

2Thanks  to  Albert  Kirk  for  making  these  tests  available. 

486 


34th  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting 


The  Allan  deviation  appears  to  confirm  that  all  generators  behave  like  flicker  FM.  The  PPL  line  is 
exactly  flat.  For  the  other  generators,  we  see  the  expected  minor  deviations  from  flatness  for  small 
r  and  an  insignificant  droop  by  DS  and  IR1  for  large  r. 


5.2  Two-Point  MSTIE 


Suppose  that  the  phase  x  (t)  of  a  clock  is  measured  at  times  to  —  t\  and  to-  For  the  purpose  of  this 
discussion,  the  mean  square  time  interval  error  (MSTIE)  of  a  clock  after  a  delay  r  is  defined  by 


MSTIE  (r,n)  =  E 


x  (to  +  r) 


x  (t0)  + 


T 

—x(to 

Tl 


(8) 


which  is  the  mean  square  error  of  linear  extrapolation  from  the  two  phase  measurements.  This 
is  independent  of  to  for  all  processes  x  ( t )  with  stationary  second  increments,  including  all  models 
and  generators  discussed  here  except  IR.  Although  the  method  of  phase  calibration  is  crude,  this 
measure  serves  the  purpose  of  showing  how  the  variance  of  the  long-term  phase  deviations  grows 
as  we  go  farther  and  farther  from  a  fixed  calibration  interval.  (See  [15]  for  a  discussion  of  more 
sophisticated  calibrations.) 

Figure  4  plots  theoretical  and  measured  values  (averaged  over  10,000  trials)  of  MSTIE(r,  10)  /r2 
against  r  for  the  flicker  FM  generators.  For  IR1,  to  =  1;  for  IR2,  to  =  1,025.  The  MSTIE  for  the 
quartz  oscillators  (normalized  as  before)  was  measured  by  averaging  the  squared  extrapolation  error 
over  all  available  to,  just  as  one  does  when  estimating  Allan  variance.  All  the  curves  except  IR1 
show  the  same  asymptotic  [1  +  In  (t/ti)]  behavior  that  has  been  calculated  for  the  PPL  model 
[?],  with  a  tiny  droop  at  the  largest  r  for  the  DS  generator.  The  IR2  curve  cannot  be  distinguished 
from  the  FD  curve,  but  the  IR1  curve  falls  off  significantly  from  the  others  as  t  increases. 


6  CONCLUSIONS 

According  to  the  tests  that  we  have  applied,  the  DS,  FD,  and  PPL  generators  are  good  fits  to  the 
phase  noise  of  a  quartz  oscillator  in  its  flicker  FM  range,  whereas  the  IR  generator  has  significantly 
smaller  long-term  mean  square  phase  deviations.  Although  the  Allan  deviations  of  the  Oscilloquartz 
and  C-MAC  oscillators  (Figure  3)  rise  as  r  increases  beyond  64  s,  this  rise  is  not  enough  to  explain 
why  the  oscillator  MSTIE  points  (Figure  4)  tend  to  line  up  with  the  curves  for  DS,  FD,  PPL,  and 
IR2  (the  “second-half”  modification  of  IR),  but  not  with  the  IR1  curve. 

The  deficiency  of  the  IR  generator  is  exposed  by  the  two-point  MSTIE,  but  not  by  the  Allan 
deviation.  Flatness  of  Allan  deviation  is  a  necessary  criterion  for  a  flicker  FM  generator,  but  it  is 
by  no  means  sufficient. 

The  FD  and  PPL  generators  are  exact  simulators  of  their  nonstationary  target  models.  Even 
though  their  outputs  live  on  a  finite  time  interval,  they  contain  contributions  from  arbitrarily  low 
Fourier  frequencies.  The  DS  generator  output,  though  it  is  only  a  stationary  process,  still  behaves 
substantially  like  a  nonstationary  /~3  process  on  its  design  interval.  We  may  regard  the  IR  gener- 
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ator  as  an  inaccurate  simulator  of  the  FD(3/2)  model.  Fortunately,  one  can  still  get  an  accurate 
FD(3/2)  simulation  from  IR  by  generating  twice  as  many  points  as  one  needs  and  keeping  only 
the  second  half.  In  this  case,  one  has  to  use  FFTs  of  size  AN  to  generate  N  points.  (One  can  also 
improve  the  accuracy  of  the  DS  generator  by  using  an  FFT  size  of  4IV  instead  of  2IV.) 

The  IR  generator  is  deficient  because  it  neglects  the  past  of  the  FD(3/2)  process;  the  same  is  true 
for  the  fractional  integration  models  in  [2]  and  the  filter  cascade  generator  with  zero  initial  state 
[9].  The  IR  output  (6)  is  just  the  first  term  of  the  right  side  of  the  FD(3/2)  formula  (2),  whose 
other  terms  represent  the  effect  of  the  random  shocks  Uj  from  the  past,  j  <  0,  on  the  future  value 
xn.  It  is  one  thing  to  tie  the  present  to  zero,  as  we  often  do;  it  is  another  thing  to  neglect  the  past 
of  these  long-nrenrory  processes. 

Finally,  despite  being  able  to  produce  models  and  generators  that  have  the  right  sort  of  behavior, 
the  author  does  not  pretend  to  understand  the  physical  mechanisms  behind  flicker  noise. 
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Fig.  1.  Phase  samples  from  flicker  FM  generators  and 
quartz  oscillators.  Generators:  DS  =  discrete  spectrum, 
IR  =  impulse  response,  FD  =  fractional  difference, 

PPL  =  pure  power  law.  Oscillators  (normalized): 

OQ  =  Oscilloquartz,  C-MAC. 


frequency  f 


Fig.  2.  Spectra  of  the  DS  generator  (N  =  32),  the 
FD(3/2)  model,  and  the  sampled  PPL  model. 

The  spectral  densities  are  multiplied  by  (2nff. 


Fig.  3.  Allan  deviation  of  flicker  FM  generators 
and  quartz  oscillators  (normalized).  IR1  =  first 
half  of  IR  output  with  N  =  2048,  IR2  =  second  half. 


Fig.  4.  Mean  square  time  interval  error,  divided  by  x2, 
for  delay  x  and  calibration  period  Xj  =  10.  All  the 
flicker  FM  generators  except  IR1  behave  like  the 
normalized  oscillators  for  large  X. 
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Questions  and  Answers 

DAVE  HOWE  (National  Institute  of  Standards  and  Technology):  Chuck,  what  is  the 
rationale  for  reflecting  the  autocorrelation  function  and  then  extending  it,  the  embedded  circulant 
approach?  Because,  it  looks  like  what  you  are  assuming  is  that  the  data  that  you  have  is  a  piece 
of  something  that  is  periodic. 

CHARLES  GREENHALL:  You  are  only  trying  to  simulate  N  plus  1  points  exactly.  What  it 
says  is  that  there  is  a  two-end  periodic  stationary  process  whose  joint  distribution  is  exactly  equal 
to  N  plus  1  successive  values  of  the  stationary  process  you  are  trying  to  simulate. 

HOWE:  My  question  is  how  do  you  know  a  priori  that  this  is  a  noisy  quasi-periodic 
phenomenon?  Why  would  you  do  the  reflection  in  a  physical  sense? 

GREENHALL:  It  is  hard  to  explain.  If  you  go  through  the  mathematics,  you  just  find  that  it 
works.  But  it  doesn’t  always  work.  It  fails  if  the  phony  spectrum  you  get  out  of  that  periodic 
autocorrelation  has  some  negative  values.  But  it  can  be  proven  that,  for  the  FD  minus  1/2  and  for 
the  other  model  that  I  worked  with  in  the  paper,  it  always  does  work. 

HOWE:  Because  ordinarily  you  would  extend  the  autocovariance  so  that  it  didn’t  appear 
periodic,  but  rather,  for  example,  you  have  truncated  data  -  you  might  assume  that  it  would 
exponentially  decay. 

GREENHALL:  We  are  not  making  any  assumptions.  We  are  just  trying  to  get  an  exact  - 

HOWE:  Right.  I  understand  that  you  are  trying  to  generate  flicker.  Second  comment:  As  far  as 
the  Allan  deviation  as  a  test  for  flicker,  you  commented  that  it  was  not  a  good  test  for  flicker. 

GREENHALL:  Well,  it  is  a  necessary  condition,  but  not  sufficient. 

HOWE:  Right.  But  in  long  term,  you  know  that  the  distribution  isn’t  symmetric.  And  is  the 
fact  that  the  last  point  drops  off  a  result  of  lack  of  number  of  degrees  of  freedom  at  the  last  point? 

GREENHALL:  Are  you  talking  about  IR1? 

HOWE:  Yes. 

GREENHALL:  No,  it  has  nothing  to  do  with  that.  That  is  just  expected  values.  That  is  theory. 
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